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Chapter 1 


1. Introduction 

Silicon is at the heart of the microelectronics revolution. Its dominance is intimately 
tied to the (i) superior materials and processing properties, and (ii) development of vast base 
of technology [1], Silicon, however, exhibits an indirect band gap and is a poor light 
emitter. For this reason, it has not found much acceptance in optical applications [1]. The 
importance of developing a technology that allows integration of optical and electronic 
devices on a silicon wafer has long been recognized [1]. Such an advance is likely to have a 
significant impact on display, communications, computer and a host of related areas. Today, 
light-emitting devices are fabricated almost exclusively from semiconductors such as gallium 
arsenide (GaAs), indium phosphide (InP), etc. which display a direct-band gap and 
relatively higher optical efficiency. Alternatively, efforts are made either to improve the 
efficiency of silicon itself or to develop silicon-compatible material with high optical 
efficiency. Considerable research has been carried out in this direction using techniques that 
range from the engineering of superlattices and quantum wells (composed of silicon, 
germanium and carbon) to the doping of silicon with optically efficient rare earths, such as 
erbium [1]. Nevertheless, the progress made is too little to improve the optical efficiency 
significantly. 

In 1990, Canham [2] observed room-temperature photoluminescence in the visible 
range from the porous silicon (PS) formed on a silicon wafer. This property of porous 
silicon was of interest for several reasons. Firstly, the emitting light had energy well above 
the band gap of bulk silicon. Secondly, the energy (or colour) could be tuned throughout 
the visible spectrum by changing the preparation conditions of porous silicon-an important 
consideration for display devices involving red, green and blue colours. Finally, the quantum 
efficiency of PS was found to be ~ 10%, i.e., comparable to that of direct-band gap 

compound semiconductors. 

1.1 Porous Silicon 

It is simply a network of nanometer-sized silicon crystallites surrounded by voids. It 
is usually prepared by electrochemical etching/anodization of silicon wafers and 
characterized by a structure consisting of a large number of micro-(diameter < 2nm), meso- 
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(2 - 50 nm) and macro-( > 50 nm) pores. Although, interest in the light-emitting properties 
of porous silicon is a rather recent development, it was discovered in 1956 during 
electropolishing ot silicon [3], A brittle, spongy, thread like structure was obtained 
exhibiting a variety of pore geometries within an otherwise crystalline material. The 
morphologies varied from a network ot fine interconnected pores to linear and wide pipe- 
like pores. Also, the pore size and shape were found to depend upon the nature of doping 
(n- or p- type and concentration), the anodization current density, electrolyte and its 
concentrations, etc. The properties of PS layers were studied extensively because of their 
potential use in electronic isolation technologies, such as silicon-on-insulator (SOI) and fully 
isolated porous silicon (FIPOS). 

PS can be prepared by several methods. The most commonly used being the 
electrochemical etching of crystalline silicon (c-Si) in a dilute solution (~ 25%) of HF and 
ethanol [1,2]. Other methods include spark erosion [4,5] and light induced etching of c-Si in 
HF [6,7], Generally, p- or n- type c-Si wafer is made an anode and etched selectively by 
passing a constant current in an electrolytic cell containing HF solution. It is indeed 
necessary to illuminate n-type substrate during the etching process for yielding porous 
silicon [8], But, a p-type wafer gives rise to PS even if etching is performed in dark. A 
variety of PS layers with a host of pore morphologies and controllable thickness can now be 
prepared. Also, successful attempts have been made to even prepare free - standing films of 
PS [9], The discovery of photoluminescence in PS at room temperature in red-yellow 
region marked the opening of an altogether new field of research because of its 
technological potential. Consequently, research work focussed around explaining the 
phenomenon of visible luminescence itself. Bulk crystalline silicon (or c-Si) has a relatively 
small but indirect energy band-gap (~ 1.12 eV), making interband radiative recombination 
process inefficient and difficult. At room temperature, c-Si shows a weak 
photoluminescence (PL) in the ER. region at about 1. 12 eV with a full width at half maximum 
(FWHM) varying with excitation conditions, but is typically ~ 0.1 eV [10,11], In contrast, 
PS shows a strong PL in the visible region at about 1. 4-2.0 eV with a comparatively narrow 
FWHM; the exact peak position, however, depends upon the preparation conditions and 
post-etching chemical treatments [2,12-14]. 

The structural aspects of PS are viewed in terms of its crystallinity which deteriorates 
with increasing level of porosity and eventually leads to an amorphous phase [15], Canham 
[2] attributed the luminescence phenomenon to widening of band gap due to quantum 



Introduction 


confinement of electrons/holes in the nanosized columnar structures. This arguments was 
based on simple Heisenberg principle and involved downward and upward shif$a of the 
valence and conduction-band edges, respectively (Fig. 1.1). Further support to the above 
explanation was provided by Raman spectra of PS, which, invariably contained a low- 
frequency side band corresponding to a shorter coherence length of optical phonons in the 
direction perpendicular to the pore axis [16], The transmission electron microscopic studies 
also reveal the presence of nano-sized crystallites in PS [17], Moreover, quantum 
confinement model has successfully explained the visible luminescence phenomenon in 
various systems containing silicon nanocrystals but produced by techniques other than 
etching [18-22], Nevertheless, Brandth et.al. [16] put forward an alternative explanation, 
now known as Siloxene model. Accordingly, visible luminescence in PS is not an intrinsic 
properties of nanostructures of c-Si, but, results from the compounds (siloxene, SUChH^; 
and its derivatives) formed by the chemical reactions of silicon with hydrogen and oxygen. 
The basis for this proposition lies in a striking resemblance found in the luminescent 
behaviour (e.g., life-time, intensity, position, shift, etc.) of PS and vibrational properties of 
Kautsky and Wohler compounds derived from siloxene [16], For example, (i) etching of PS 
with an electrolyte containing of mixture of HF, HC1, and ethanol results in a luminescence 
band centred at 690 nm. This agrees well with the fluorescence observed in siloxene after 
annealing at 673 K, the larger peak width here is due to rather uncontrolled substitution of 
hydrogen by OH groups during the treatment, (ii) the red luminescence in PS found after 
etching in a solution of HF and ethanol has a maximum at 760 nm, which is in quantitatively 
agreement with the fluorescence peak of siloxene after annealing at 673 K. Infra-red (ER) 
and Raman spectra further corroborated to the siloxene model [16]. In yet another model, 
etching process is believed to form Si-hydride on the large surface area of PS and excitations 
in them are responsible for the photoluminescence [23,24], Moreover, the distribution of 
hydrides determines the peak position while total surface area of PS layers influence the 
intensity of the PL spectra. The experimental observations that led to this model include red 
shift and gradual disappearance of PL on annealing at temperatures in the range of 503-663 
K, blue shift in PL peaks due to pore widening in PS after chemical etching with HF [23], 
etc. It can therefore be inferred that photoluminescence (PL) phenomenon in porous silicon 
is still not fully understood. 

Apart from visible PL, electroluminescence (EL) was observed in PS though with 
disappointingly poor efficiency (~ 0.01%) [25-27], This initiated studies on PS with regard 




1 . 1 Schematic diagram of a semiconductor quantum well. 



to carrier transport properties too. Also, metal contacts on PS layers are found to act as 
light emitting diodes [25], It is worth mentioning here that PS in general has a complex 
microstructure : quantum wire columnar structure, oxidised (or otherwise passivated) 
partially connected spheres of silicon with tissue like morphology, amorphous with small 
embedded Si nanocrystals, etc. Its exact nature, however, varies with the substrate 
resistivity, porosity level, and thickness. 

1.2 Objective of present work 

Ever since visible PL was discovered in porous silicon [2], enormous studies have 
been undertaken using different characterization techniques for obtaining complimentary 
information [ 28-3 1], However, thermal stability aspects and the nature of chemical species 
present in PS have not been investigated in depth [32], Since the mode of preparation of PS 
involves electrochemical bath containing HF and alcohol, many different chemical species 
are likely to be present. Further, PS has a large surface-to-volume ratio (~ 4 x 10 8 m 2 /m 3 ) 
and very high density of dangling bonds at the surface [33 ]. These dangling bonds may act 
as traps or non-radiative recombination centres for the photo-generated electron-hole pairs. 
There exists several reports on optical and transport properties of the PS to suggest that the 
processes at the surface play a dominant role in determining its behaviour [34-36], It is 
important therefore to have a reliable and quantitative method of identifying surface 
desorbed species in PS for evaluation of its thermal stability. In the present work, an 
attempt is made to prepare porous silicon by usual electrochemical process and monitor its 
effusion characteristics on heating in vacuum ~ 1CT 6 torr by thermally stimulated effusion 
(TSE) and FTIR spectroscopy. In addition, complimentary techniques consisting of Raman 
and e.p.r. spectroscopies and X-rays diffraction have been utilized ^gather information with 
regard to microcrystallinity, dangling bonds/spin density, phase identification with 
orientation disposition and grain size, etc. for arriving at the better . insight of the nature of 
porous silicon. 



Chapter 2 


2. Experimental details 

This chapter contains details of the preparation of porous silicon and 
characterization techniques used. Sample preparation is a long process and involves 
several steps, viz., cleaning of silicon wafers, metal deposition for back contact, 
electrochemical etching under controlled conditions, drying, etc. For characterization of 
porous silicon, Raman spectroscopy. X-ray diffraction, electron para magnetic resonance 
(EPR) spectroscopy, thermal gas effusion and fourier transform infrared (FTIR) 
spectroscopy are employed. 

2.1 Formation of porous silicon 

2.1.1 Substrate Preparation 

The substrate used was a p-type (boron doped) polished c-Si wafer of resistivity p 
= 0.5-10:Q-cm and orientation <100>. It was cleaned chemically in a multistep process. 
First the wafer was boiled in trichloroethylene (TCE) to degrease and rinsed with ethanol 
to remove the residue. The wafer was further cleaned in an ultrasonic cleaner using 
acetone. The oxide layer that usually exists on Si-wafer was removed by etching using a 
dilute solution (20%) of HF in ethanol for about 30-40 s. It was then rised with deionised 
water. The resulting surface was generally hydrophobic in nature. For back contact, 
aluminium was evaporated thermally in vacuum (~ 10' 5 torr) and allowed to deposit on the 
unpolished (or back) surface of c-Si wafer. The evaporation source was a tungsten 
filament (i.e., a basket) and held underneath and facing the back surface of c-Si. To ensure 
high quality of dWmiuw film, evaporation was carried out under vacuum ~ 10' 5 torr or 
better. Due to the oxide surface layer invariably formed during the transfer of c-Si into the 
vacuum chamber even, there develops a barrier across the silicon-aluminium interface. 
This makes the contact non-ohmic in nature. To overcome this difficulty, the aluminium 
coated silicon was annealed at ~ 550°C for about half an hour. At this temperature, 
interatomic diffusion occurs at the interface leading to alloy formation. Fig 2.1 shows the 
schematic diagram of the annealing set-up. It consists of a hollow cylindrical furnace, a 
quartz tube connected to a pumping system and arrangement of measuring the temperature 
inside through a thermocouple. Aluminium back coated c-Si wafers were placed inside the 
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Fig. 2.1 : Schematic diagram of annealing set-up 


attained, the quartz tube was inserted into the furnace already maintained at the set 
temperature. After annealing for 0.5 h, the furnace was turned off and sample was left to 
cool down as such to room temperature. Finally, the pumping was shut down and the 
sample taken out for electrochemical etching. Fig 2.2 gives the flow diagram of substrate 
preparation steps. 

2.1.2 Electrochemical etching 

Fig 2.3 shows a schematic diagram of the electrochemical cell used for etching of 
c-Si. It consists of a teflon cylinder with a coaxial circular slot engraved for holding a 
mask to enable selective exposure of c-Si wafer. The mask is just a washer made of teflon. 
An O-ring is put inside the groove (around the central hole) to seal and prevent acid 
leakage. For etching, c-Si wafer is placed in the central slot on the mask with the face to 
be etched in contact. A brass disc fitted coaxially with a rod is polished well using a fine 
grade emery paper and is placed touching the aluminium coated c-Si for making electrical 
connection. Threads are made on the walls of teflon cylinder to fit a nylon nut having a 
through hole to accommodate the brass electrode. By tightening the nut, optimum 
pressure can be applied on the brass electrode to obtain a good electrical contact with the 
c-Si wafer. The entire assembly is placed into a nylon seat with the polished c-Si surface 
facing upwards. The nylon seat is also a hollow cylinder with a slot of diameter matching 
with the outer diameter of the teflon cell. The lower end of the nylon cylinder is provided 
with a notch for leads to complete the electrical circuit. 


Table 2. 1 . Preparation conditions of porous silicon (PS). 


Sample 

code 

Resistivity of 
Si-substrate 
p(D-cm) 

Etching Current 
(mA/cm 2 ) 

Etching time 
(minutes) 

PSL1 

0.5-1. 0 

10.5 

120 

PSL2 

0.5-1. 0 

10.5 

90 

PSL6 

0. 5-1.0 

5.0 

120 

PSL11 

0.7-1. 5 

5.0 

15 

PSL13 

0.3-0. 6 

10.0 

30 
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Fig. 2. 2 :Flow diagram of the substrate prepration steps. 
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Fig2.3:A set-up for electrochemical etching of silicon. 


The 1:1 mixture (by volume) of 48% HF and 99.9% C2H5OH is poured into the 
teflon cell. The whole set-up is mounted on a stand. A platinum wire immersed into the 
electrolyte serves as a cathode. A constant current is passed through the cell for 
anodization to occur at the exposed surface of the c-Si wafer. While Kiethley 220 
programmable current source is used as the power supply, a Keithley 196 DMM monitored 
the voltage across the cell. The whole set-up is housed in a clean-bench to not only 
provide a dust free environment, but also, safeguard against the hazards of HF fumes. 
Initially, a high current (raising the density to 100 mA/cm 2 ) is passed for about 10 s to 
electropolish and make the surface uniform. The current source is programmed before 
resuming etching. As the process begins, bubbles emerge at the substrate surface. 
According to Smith and Collins [37], SH2 gas evolves during bubbling and pore formation 
occurs in c-Si top layers. After etching, electrolyte is removed and both the c-Si wafer and 
the cell rinsed with de-ionized water for 10 minutes and dried. For a given concentration 
of HF, thickness of porous layer formed is determined by etching time itself A beam of 
UV light (from 125 W mercury lamp) is then impinged onto the etched area to examine 
luminescence, if present. The resulting porous silicon samples indeed exhibited the typical 
red or yellowish red luminescence. The samples are stored in a desiccator in dark. 

2.2 Characterization methods 
2.2. 1 Raman scattering 

Raman effect basically refers to the emergence of radiation of 
wavelength/ffequency other than the incident beam at some angles on interacting with 
certain molecules. The process involves inelastic scattering and is very sensitive to the 
vibrational modes of molecules/atoms. Also, the extent of change (or shift) in frequency 
depends upon the chemical structure of the molecules responsible for the scattering. 
Raman spectrum is commonly plotted in terms of wavenumber shifts with respect to the 
source (incident radiation) line. The shifts towards lower and higher energies are termed 
as stokes and antistokes; the intensity of antistokes line is appreciably smaller than the 
corresponding stokes lines. As discrete energy shifts are observed in the Raman spectra, 
they indicate the presence of quantized energy transitions. The shift is understood by 
assuming that the incident radiation on interaction with the sample causes polarization and 
energy retained momentarily (i.e., for 10' 15 - 10" 14 s) in the virtual state without involving 
transition to a higher electronic energy level. Fig 2.4 shows interaction of radiation with 
zero and first vibrational levels of ground electronic state. After retention, the species 



return to their ground state with release of radiation in all direction with exactly the same 
energy as the source (shown by dashed arrows in Fig 2.4). However, the species may relax 
to the first vibrational level of the ground state (broad arrows in Fig 2.4). The frequency 
of light emitted in this transition will be less by a quantized amount (AE), i.e., difference in 
energies of ground and first vibrational levels. Thus, 

AE= h (Av) 

where Av is the frequency shift observed in the Raman spectra and h is Planck’s constant. 
On the other hand, when the molecule already in the first vibrational level goes to the 
virtual state by interacting with the incident radiation and return to the ground state will 
release energy more by an amount AE. Since the probability of molecules being in first 
vibrational state is low than in ground state at room temperature, antistokes radiation is 
weaker than stokes. 

It may be mentioned that both the Raman and the infrared spectra depend upon the 
same vibrational mode but arise from processes that are mechanistically different. While 
infrared absorption requires a change in dipole or charge distribution of the molecule and a 
direct transition, Raman scattering involves a momentary polarization of the molecule that 
ceases upon relaxation and emission. The effectiveness of Raman scattering thus depends 
upon the case with which electrons of the bond can be distorted from the normal positions 
(i.e., polarizability of the bonds); polarizability decreases with increase in electron density 
and bond strength but decrease in bond length. The Raman shift, in fact, requires a change 
in polarizability rather than change in dipole moment - associated with the vibrational mode 
of the molecule. For example, the molecule such as nitrogen, chlorine or hydrogen has no 
dipole moment and so not active towards infrared absorption. But, the polarizability of the 
bond between the two atoms varies periodically in phase with the stretching vibrations 
reaching a maximum and minimum at the extreme and closest separations, respectively. A 
Raman shift corresponding in frequency to that of the vibrational mode thus results. 
Raman spectrum gives a spectroscopic signature of structure, symmetry and dynamical 
characteristics of molecules. 

Porous silicon samples were examined in a SPEX-1877 Triplimate Raman 
spectrometer. For this, 20 mW light of 514 nm wavelength emanating from an Ar" laser 
was focussed onto the sample with a spot size of 50 itm and back scattered rays were 
collected by a lens, passed through a slit to a triple monochrometer and finally to a 
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multichannel detecting system, consisting of charged coupled device (CCD) and photon 
counter. The CCD is cooled to liquid nitrogen temperature. The arrangement leads to a 
high sensitivity (quantum efficiency > 60%) with an extremely low background and noise 
level, thus enabling long accumulation time. The spectrum data is fed to a personal 
computer (PC) for storage, display and further analysis. 

2.2.2 X-ray diffraction 

Rich Siefert ISO - Debyeflex 2002 X-ray diffractometer with CuKot radiation 
X = 1.5418X was employed to record the diffraction patterns and identify the phase (s) 
present in the porous silicon samples and changes occurring on annealing at 250°, 450° and 
650°C. The calibration of the system was carried out with a standard silicon sample. 
Further, the average size (t) of crystallites present in the sample is evaluated from the X- 
ray line broading [38], Accordingly, the average crystallite size is given by the Schemer 
relation 


t = — — with B = B -B s 

BCosO 

where X is the wavelength of X-ray beam, B is the FWHM of diffraction peak (measured in 
radians) at a particular Bragg angle 8 (corresponding to 20 in the diffractogram); B m and 
B s are FWHM of diffracted peaks obtained with the sample and the standard (c-Si, in the 
present case), respectively. The standard sample having a grain size > 0.2pm is used to 
estimate the instrumental broadening and making correction in the experimental data. 

2.2.3 EPR spectroscopy 

Electron paramagnetic resonance (EPR) spectroscopy is based upon the absorption 
of microwave radiation by solids when placed in a strong magnetic field. Species that 
contain unpaired electrons are responsible for the absorption and therefore permit their 
detection. The electron has two energy levels, lower with spin quantum number m = -1/2 
and higher with m = +1/2. These energy levels split under influence of magnetic field (Ho), 
such that their difference (AE) is given by 


AE = g (3 n Ho 



where ‘g’ is splitting factor and p N is Bohr magneton (=9.27 x 1CT 21 erg G' 1 ). The value of 
‘g for free electron is 2.0023 but varies with environment, e.g., is a few % of 2.0023 for 
unpaired electron in a molecule or ion. The value of AE is 6.3 x 10' 24 J or 0.03 meV in a 
magnetic field of 3400 gauss. This energy difference corresponds to frequency (e.g., 9500 
MHz) that lies in the microwave region. The microwave radiation of frequency (say 9.5 
GHz) emanating from source falls onto the sample held in a small quartz tube positioned 
between the pole of the permanent magnet. Helmhotz coils provide a means for varying the 
magnetic field over the small range (~4 Gauss, 100 KHz) in which the resonance occurs. 
The absorption of radiation can then be measured as a function of magnetic field. 
Generally, EPR spectra are recorded in derivative form of the Intensity vs. magnetic field 
data to enhance sensitivity and resolution. A Varian model E-109 spectrometer has been 
used in X-band to record the derivative microwave absorption spectra of virgin porous 
silicon samples (PSL 2 and PSL 6) and after annealing them upto 250°, 450° and 650°C. A 
standard sample of ococ'-di phenyl-p-picryl hydrasyl (DPPH) is employed to determine the 
‘g’ factor of porous silicon using 

^n(DPPH) 

(f — CT Y 

* PS ~ 0 > DPPH ^ rr 

** o(PS) 

where gm>m = 2.0036 at room temperature, and H 0 (dpph) and H 0 (ps> are resonance magnetic 
field found in the EPR spectra of DPPH and PS, respectively. The signal parameters (viz., 
peak to peak width and height) are measured from the EPR spectra for the estimation of 
spin density as described later in chapter 3. 

2.2.4 Thermal gas effusion 

Thermal gas effusion is employed to identify chemical species released on heating 
of porous silicon (PS) in vacuum. The set-up used for this purpose is shown schematically 
in Fig. 2.5. The porous silicon (PS) sample is placed in a quartz tube (inner diameter 1.44 
cm and length - 18 cm), connected to a quadruple gas analyzer (QGA) and a turbo 
molecular pump (TMP). For temperature measurement, a chromel-alumel 
thermocouple is inserted into another quartz tube of inner diameter of 3 mm and placed by 
the side of the porous silicon (Fig. 2.5). Vacuum ~ 10‘ 6 torr is created inside the quartz 
tube. The temperature is then raised to a particular level by switching on the furnace. The 




Fig. 2.5: Schematic diagram of Thermal Effusion Set-up 
DPM: Digital Panel meter DPH: Differential 






gasa evolved are analysed with a QGA by converting them into respective ions and 
measuring their m/e ratios. 

A quadrupole gas analyzer (QGA) contains four short parallel metal rods arranged 
symmetrically around the ion beam (Fig. 2.6). The opposite rods are connected together, 
one pair is attached to the positive while the other to the negative of a d.c. power source. 
In addition, a radio frequency a.c. potential is applied to both the pairs. Neither field acts 
to accelerate the positive ions. The combined field, however, causes species to oscillate 
about their central axis of travel; only those with a certain mass-to-chayrge ratio can pass 
through the array without being removed by collision with one of the rods. Mass scanning 
is achieved by varying the frequency of the a.c. supply or changing the potential of the two 
sources while keeping their (frequency/potential) ratio constant. The ions which pass 
through are collected by an electrode, well shielded from stray ions, and current thus 
produced is amplified electronically and recorded. Ion current usually vary over a range of 
10’ 17 to 10'%. QGA works in two modes : (i) the partial pressure of gases versus m/e ratio 
and (ii) partial pressure of gases (12 in number) as a function of time. In the present work, 
QGA is used in the first mode for various ions listed in Table 2.2- A personal computer 
(PC) was interfaced with QGA for monitoring the partial pressures as a function of 
temperature and storing the data. 


Table 2.2 The ions traced in thermal effusion of porous silicon. 


Serial no. 

(m/e) of ions 

Ions 

I. 

1 

H + 

2. 

2 

h 2 + 

3. 

16 

0 + 

4. 

29 

SiHT 

5. 

30 

SiH 2 + 

6. 

31 

SiH 3 + 

7. 

32 

o 2 + 

8. 

47 

SiF + 

9. 

48 

SiFHf 

10. 

49 

SiFH 2 + 

11. 

68 

SiF 2 + 

12. 

87 

SiF 3 + 





FIGURE 



2.-6 \ quadrupole mass spectrometer 
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2.2.5 FTIR spectroscopy 

The fourier transform infrared (FTIR) spectrometer is shown schematically in Fig. 
2.7. The radiation from a broad-band light source is divided into two parts by a beam 
splitter. After reflection from a fixed and a movable mirror, respectively, the two partial 
waves are combined by the same beam splitter and focussed on a detector. They interfere 
there and yield the total intensity. For different positions of the movable mirror the two 
partial waves get different phase shifts with respect to each other. Therefore, on the 
detector the radiation field is superimposed with a time-delayed copy of itself. Hence, the 
detector signal resulting while the mirror moves is basically is the autocorrelation function 
of the radiation field (or interferogram). The fourier transform of this function gives the 
FTIR spectrum in the frequency domain. The spectral resolution in FTIR is determined by 
the inverse of the total path of the movable mirror. The light from all spectral intervals 
contributes to the detector signal and no mechanical slits are required here. These features 
reduce the time needed to record a spectrum considerably. FTIR spectra were recorded in 
the laboratory atmosphere using a 1760 - X Perkin - Elmer spectrophotometer in the 
wavenumber range of 500 - 3600 cm' 1 with a resolution of 0.1 cm' 1 . As aluminium itself 
absorbs infrared radiation, measurements were carried out after removing that from the 
back of the porous silicon. FTIR spectra were recorded for virgin samples (PSL2, PSL6 
and PSL13) and again after annealing them upto 250° and 450°C. 

2.3 Peakfit Analysis 

For analysing the Raman spectra and thermal gas effusion data, a software package 
called “peakfit analysis” was employed. It is a menu driven package and therefore easy to 
use. In this, the data is first plotted, and then subjected to base-line correction and curve 
fitting with distributions such as Gaussian, Lorentzian, Voit, etc. Subsequently, by 
adjusting the width and the height of peaks, iterations performed until convergence is 
attained. The fitted curve then provides information about peak position, full width at 
half maxima (FWHM), percentage area under various peaks, etc. as discussed in chapter 3. 
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Fig. 2.7 Scheme of a Fourier transform spectrometer 
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Chapter 3 


3. Results and Discussion 

In this chapter, we present the characterization results of porous silicon involving 
Raman, X-ray diffraction, electron paramagnetic resonance (EPR), thermal gas effusion 
and fourier transform infrared (FTIR) spectroscopic techniques. Also, the effect of 
annealing on the nature of PS sample is described using X-ray, EPR and FTIR data. 

3.1 Raman Spectroscopy 

Fig. 3.1 shows a typical Raman spectrum of p-type c-Si together with the curve 
joining experimental data points. Its analysis carried out using a peakfit software package, 
as mentioned in chapter 3, gives deconvoluted plot depicted by dashed curve (not visible 

since superimposed on the experimental data curve of Fig 3.1). The central maxima occurs 

a 

atJRaman shift of 520. 1 cm' 1 with a FWHM of 7.0 cm' 1 . It can be noticed that the peak is 
symmetric in nature. Also, the values obtained match well with the reported Raman data 
of c-Si wafers, i.e., the peak maxima and FWHM lying in the wavenumber ranges 
520-522 cm' 1 and 3-8 cm' 1 , respectively [31,39-41], According to Iqbal and Vepek [30], 

silicon films exhibit Raman spectra with peak maxima at progressively lower frequencies 

in 

and larger FWHM besides increase^asymmetry (i.e., tail towards low frequency side) with 
decrease in grain size beiow 200 nm. Raman spectrum of an a-Si contains usually a very 
broad and weak peak near 480 cm' 1 and is attributed to the vibrational density of states 

[ 30 ]. 

Raman spectra of virgin porous silicon samples PSL2, PSL6 and PSL13 recorded 
from two surface regions (say ‘a’ and ‘b’) of diameter 50 pirn each are depicted in Fig 3.2- 
3.4. While the experimental profiles are drawn by continuous curves, deconvoluted plots 
obtained for each are shown by three dashed Gaussian distribution functions. Table 3.1 
gives summary of three peak maximas and their FWHM ranges for various Raman spectra 
(Fig 3. 1-3.4). As can be noticed, the first peak (No.l) maxima lie within the wavenumber 
range 519-521 cm' 1 and therefore correspond to c-Si. However, widths of signal peaks 
emerging from region ‘a 3 are relatively larger than from region *b\ which actually 
corresponds to more of crack portion and hence give contribution of the substrate itself. 
The larger FWHM of Raman peaks resulting from region ‘a’ can be attributed to both the 
microcrystallites present in PS and c-Si substrate. Red shifted peak (No. 2) within 
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wavenumber range 502-513 cm' 1 having a larger FWHM of 23-28 cm' 1 is due to the 
nanocrystallites (< 5nm) and therefore considered as charactersitic of porous silicon [41], 
These results indicate that porous silicon region is not homogeneous with regard to • 
distribution of crystallites and pores. A very broad peak (No. 3) having FWHM in the 
range 48-61 cm' 1 and centred around wavenumber 473-490 cm' 1 appears because of the 
presence of disordered a-Si phase in porous silicon [39-41], Asymmetry found in the 
Raman spectra towards lower wavenumbers itself suggest the existence of nanosized - Si 
crystallites in PS. It needs to be mentioned that annealing of PS makes the Raman spectra 
progressively symmetrical in nature as the crystallites grow in size in the process [42], 


Table 3.1: Peak positions and their FWHM ranges for various Raman spectra shown in 
Figs. 3. 1-3.4 of porous silicon. 


Sample (location) 

Peak 1 

Peak 2 

Peak 3 

Position 

(cm' 1 ) 

FWHM 

(cm' 1 ) 

Position 

(cm' 1 ) 

FWHM 

(cm' 1 ) 

Positio 

n 

FWHM 

(cm' 1 ) 

c-Si 


520.1 

7.0 

- 


- 

- 

PSL2 

(a) 

520.8 

15.1 

512.6 

25 

486.9 

47.7 

PSL2 

(b) 

520.3 

8.6 

508.8 

26.5 

479.1 

48.3 

PSL6 

(a) 

519.7 

10.4 

507.6 

26.5 

490.4 

60.5 

PSL6 

(b) 

520.3 

3.0 

512.2 

22.6 

- 

- 

PSL13 

(a) 

519.8 

7.1 

502.1 

27.1 

473.0 

55.3 • 

PSL13 

(b) 

520.0 

6.6 

510.5 

27.5 

479.9 

56.9 


The redshift and peak broading observed in Raman spectra are usually attributed to the 
spatial confinement of optical phonon with decrease in crystallite size [43,44], The optical 
phonons at nanosized crystallites do not correspond to plane waves and their localization 
causes relaxation in the wave vector q*0 selection rule. As a consequence, phonons of not 
only q=0 but also q>0 take part in Raman scattering process, resulting in the redshift and 
broading of Raman peak. Fig 3.5 gives the redshift of PS as a function of FWHM. 
Clearly, the presented data (as found by others as well [45]) do not reflect any trend. 
Actuallv. ohonon confinement model predicts a parabolic nature of such a plot [45], 
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Yang et. al. [45] have, however, shown that phonon confinement description do not match 
with the experimental results. Instead, effect of strain (usually ~ 10’ 3 ) is also to be 
considered to account for the experimental data. 

3.2 X-ray diffraction 

Figs 3.6 and 3.7 show the X-ray diffractograms of porous silicon samples PSL2 
and PSL6 in the 20 range of 58-78° recorded with CuKa radiation (A = 1.5418/4). All 

virgin samples exhibit a broad peak around 29 = 69.2° which corresponds to 400 reflection 
of silicon. Needless to say, silicon wafers used for preparing PS had (100) orientation. 
Therefore, one expects reflections of hOO type only and the structure factor considerations 
limit those which satisfy h=4n (n being an integer). These observations (i.e., peak broading 
and 400 reflection) indicate that crystallities belong to diamond type silicon with lattice 
parameter a=5. 4307 A . are of nanometer size and maintain the orientation of the 
original silicon wafer. The line broad-^yyiata and corresponding average crystallite size 
determined (as described in section 2.2.2) are given in table 3.2. 

After annealing at temperature at 450° and 650°C (see, e.g. Fig 3.6) of PSL2, 400 
reflection peak becomes narrower with a sharp central maxima. It appears as if these are 
two peaks, originating from (i) c-Si substrate depicted by dashed curve and (ii) porous 
silicon shown by dotted curve, superimposed. Considering narrow peak for c-Si as 
standard and broad for PS, one can make corrections in the FWHM values and estimate 
changes in the crystallite size. The values thus obtained are also given in Table 3.2. These 
clearly show an increase in crystallite size by ~33 and 52% at 450° and 650°C, respectively. 
The corresponding results with PSL6 samples reveal decrease in peak width without sharp 
central maxima found in PSL2 (Fig. 3.7). The crystallite size here increased by 25 and 
50% at 250° and 450°C, respectively (Table 3.2). The difference in behaviour is caused as 
< some portion of porous silicon region was lost in PSL2 during handling. This turned out 
be advantageous as standard became available right on the sample- itself and exactly at the 
diffraction angle 20 = 69.2°. The average crystallite size is invariably smaller in PSL6 
mainly because of the longer etching time used [46], As discussed later, effusion of 
hydrogen and other chemical species take place during annealing of PS. Consequently, 
average crystallite size should decrease. But, annealing also causes coarsening of particles. 
Since, an overall increase in crystallite size is observed, it can be inferred that process of 
coarsening predominates if effusion species rest on the surface only. Alternatively, the 
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civstallitcs must turn porous thereafter. Also, the crystallites continue to maintain 
preferred orientation (100) even after annealing atleast upto 650°C. 


I able 3.2. X-ray line broading data of porous silicon with average size of crystallites. 


Sample 

Virgin/ 

annealing 

temperature 

(°C) 

29 

(degree) 

B m 

(degree) 

B s 

(degree) 

B 

(degree) 

Average 

crystallite 

size (t) 

(nm) 


Virgin 

69.2 

5.5 

0.9 

4.6 

2.1 

PSL2 

450 

69.2 

4.4 

0.9 

3.5 

2.8 


650 

69.2 

3.9 

0.9 

3.2 

3.2 


Virgin 

69.2 

7.0 

0.9 

6.1 

1.6 

PSL6 

250 

69.2 

5.8 

0.9 

4.9 

2.0 


450 

69.2 

5.0 

0.9 

4.1 

2.4 

B m and B, are FWHM of diffraction peak of PS, and c-Si; respectively; 

B •• (B m - B„), / = A = 1.5418/d 

BCosO 



3.3 Electron paramagnetic resonance (EPR) spectroscopy 

The electron paramagnetic resonance technique is particularly suited for the studies 
of bulk, surface and interface defects in Si [47,48]. The EPR signal must come from the 
porous Si layer as it is observed in anodized samples and not in the original Si wafers. 
When the porous layer is thinned artificially by polishing, the intensity of the signal 
reduces progressively, which means that .the spin centres are associated with the porous 
layer itself. Also, the detectable signal is found not only, to be, identical in the various 
porous layers, but also, independent of the substrate resistivity, substrate orientation and 
anodization conditions [49], All these results confirm that the EPR signal is intrinsic to the 
porous Si. The following arguments are pertinent and useful for associating the EPR 
signal to Si dangling bonds. 










(a) From the previous EPR data of defects in c-Si, such as the phosphorus-vacancy 
pair [50], the five-vacancy cluster [51], the nitrogen-vacancy pair [52], and the P b centre at 
the Si/SiOj interface [53], we find that the g-tensor is either exactly or approximately 
trigonally symmetric (axial direction <11 1>) (or deviating a little, usually by ~ 3°). 
Moreover, these defects have their spin electrons localized at Si dangling bonds which are 
oriented along <11 1> directions. The slight deviation from <1 1 1 > may be caused by 
perturbation due to impurity or other lattice defects present near the dangling bond. 
Obviously, trigonal symmetric EPR signal can be identified with Si- dangling bonds in c-Si. 

(b) The principal g values obtained experimentally in respect of various defects whose 
spin electrons are localized at Si dangling bonds are listed in table 3.3. It is noticed that 
the value of gn and g L lie in ranges 2.005-2.0026 and 2.0081-2.0110, respectively. 
Watkins and Corbett et al [50] have also found theoretically gn=2.0023 (which matches 
with free electron case) and g L -2.010 for unpaired electrons associates with Si-dangling 
bonds. 


Table 3.3: Comparision of EPR g-values for spin centres consisting of an unpaired electron 

at a Si-dangling bond. 



Identity 

gn 

gi 

Refs. 

c-Si 

phosphous vacancy 

2.0005 

2.0104 

[50] 


five vacancy 

2.0020 

2.0110 

[51] 


Nitrogen-vacancy 

2.0026 

2.0089 

[52] 

Si/Si0 2 

P b center 

2.0012 

2.0081 

[54] 

Interface 


2.0013 

2.0090 

[53] 

Porous silicon 

Dangling bonds 

2.0023 

2.0086 



The effective g-values determined from EPR spectrum for the porous silicon PSL2 and 
PSL6 are found to between 2.0061 and 2.0067. In bulk c-Si , the dangling bonds may 
possess their axial directions along four equivalent <11 1> axes [50-52], while at the (111) 
Si/Si0 2 interface, they are mostly aligned along direction perpendicular to interface 
[49,53,54]. However, for (100) Si/Si0 2 interface, dangling bond direction will be incl.*™* 
to the interface normal [100], So, their contribution to the EPR signal is expected to be 
anistropic. That is perhaps the reason for our obtaining g-values in between the gn, g L 
and as effective. On the other hand, the dangling bonds in an a-Si can not have 





such an orientation disposition because of the absence of a regular or spatially coherent 
net w oik I icvious and piesent results of Raman spectroscopy have revealed an amorphous 
region near the top layer of the porous silicon with characteristic peak at 480 cm" 1 . The 
conti ihution of the amorphous state to EPR, if exists, must be isotropic in nature. 

Due to the low spin electron density in c-Si, their detection in bulk samples of a 
few mm size has not been possible. However, in the case of porous silicon, there is large 
internal sui face area and so total spin number may become adequate to give a detectable 
bPR signal I herefore, the proporties of the Pb centre can possibly be used to characterize - 
the internal surface of porous silicon. Finally, thermal annealing at 250°, 450°, and 650° 
under ultrahigh vacuum may determine the importance of hydrogen passivation of Pb 
centres in the virgin samples. Both virgin samples PSL2 and PSL6 yield EPR signals. Figs 
3.8 and 3 9 show room temperature EPR spectra of PSL2 in virgin state and after 
annealing upto 450 <I C at a rate of 9.6°C/min for 2 hours as a function of the magnetic field 
perpendicular to the sample surface (i.e., along [100], direction). It can be noticed that the 
line shapes are almost Lorentzian. Also, the line widths of the first derivative absorption 
spectra are 6.7 G and 6.0G for PSL2 and PSL6, respectively. The principal value of g 
tensor in virgin state of samples were determined as gj>si.2 = (2.0061 ± 0.0002) and gi>si.r> = 
(2.0067 ;l 0.0002). In PSL2, two peaks were observed in the EPR spectra after annealing 
at 650°C and so the corresponding g-values determined are gn = 2.0073 and g, = 2.0035. 
The effective g-values was obtained from 


g fir = + 2g, )/3 as 2.0045. The result of EPR studies are 

summarized in table 3.4. It is evident that the signal increases upon annealing in both PSL2 
and PSL6, indicating the increase in the number of dangling bonds. The spin numbers (I) 
is calculated using the formula 

I = KW 2 H 

where W is the peak width and H is peak height of signal in derivative mode and K is a 
constant. The values of K for line shapes in respect of Lorentizian and Gaussian 
distributions are 1.81 and 0.52, respectively. The standard DPPH sample is assigned a 
value of 1.1 for the constant K. Fig 3.10 shows the variation of ‘g’ values with the 
annealing temperature. It is clear that g-values change somewhat upon annealing but 
remain in between 2.0033 and 2.0067. Fig. 3.11 demostrates that the peak width 
increases continously from 6.7 to 7.5 G upon annealing of PSL2. However, the changes 




3.8 EPR spectrum ofPSL2 in virgin state and the standard DPPH sample. 




EPR spectrum of PSL2 after annealing upto 450 C and the 
standard DPPH sample. 
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are irratic for PSL6. Fig. 3.12 depicts the relative spin numbers vs. annealing temperature 
plot. Obviously, the spin number increases to 45 times for PSL2 after annealing at 450°C 
and 10 times for PSL6 at 650°C. It means that the dangling bonds are indeed created 
upon annealing. 

The thermal gas effusion data (discussed in section 3.4 ) indicate that the hydrogen 
and hydride species (SiH x ) are evolved during annealing. This may produce dangling 
bonds, resulting in an increase in the EPR signal, as observed here. These results are in 
qualitative agreement with those reported earlier [50,53], 


Table 3.4 Effect of annealing on EPR signal and related parameter in porous silicon. 


Sample 

Annealing 

temperature 

(°C) 

g-values 
( ± 0.0002) 

peak 

width 

W 

(Gauss) 

Height 

of 

signal 

(H) 

(arbitray 

unit) 

W 2 H 

Spin 

numbers 

(arbitray 

unit) 



Virgin 


6.7 

26 

1167 

1167 

1.0 

PSL2 

450 


7.5 

954 

5366 

53662 

46.0 


650 


two 

peaks 


** 

** 

“ 


Virgin 

2.0067 

6.0 

56 

2016 

2016 

1.0 


250 

2.0063 

5.0 

37 

925 

925 

0.5 

PSL6 

450 

2.0064 

6.0 

225 

8100 

8100 

4.0 


650 

2.0063 

5.5 

614 

1858 

18581 

9.5 


3.4 Thermal gas effusion 

The gas effusion from porous silicon PSL1 and PSL1 1 samples annealed at the rate 
of 9.6°C per minute under clean vacuum ~ 1C 6 torr was monitored by instant partial 
pressure at various temperatures and setting m/e numbers for H", H 2 + , 0 + , SilT, SiH 2 + , 
SiH 3 + , 0 2 + , SiF + , SiFH + , SiFH 2 + , SiF 2 + and SiF 3 + in a quodrupole gas analyzer (QGA). 
Figs. 3.13 and 3.14 show the gas effusion spectrum as a function of temperature for 
samples PSL1 and PSL11, respectively. The partial pressure of gas species is invariably 
larger in PSL1 presumably due to thicker layer of porous silicon as etching performed for 
longer time (120 m) against 15 m for PSL1 1. The assignments of various m./e numbers is 
based on the likely evolution of chemical species and their cracking patterns. For example. 
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p,q . 3.13 Thermal gas effusion spectrum of PSL1 as a function of temperature: 
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3.14 Thermal gas effusion spectrum of PSL1 1 as a function of temperature. 




m/e=2 was assigned to the molecular hydrogen and not to He 1 * as the latter is not 

expected to be present in any significant amount in porous silicon sample because of the 
use of HF solution for etching. Similary, m/e=29, 30 and 31 are attributed to SiH\ SiH 2 + 
and SiH?', respectively. Higher m/e assigned to fluorides of silicon. 

The individual spectrum was analyzed using a peakfit software package. Fig. 3.15 
shows the case of molecular hydrogen evolution from PSL1. It contains two peaks central 
around 376° and 476°C with FWHM of 38.6° and 139.7°C, respectively. Figs. 3.16 - 3.18 
depict the spectra and their fitted curves for SiH, SiH 2 and SiH 3 species. It is interesting to 
note that all of them have central maxima within 330° - 350°C. These results suggest that 
SiH species are evolved nearly at the same temperature and presumably lie on the surface 
of the crystallites. The hydrogen trapped in pores and present interestitially in silicon 

crystallites gets released around 376° and 476°C. 

Assuming that the kinetics of gas desorption processes follow a rate equation [51] 

— = ,4(1 - N) n exp (- EJRT) 0 ■ 0 

dt 

where N is the fraction of gas evolved by overcoming the free energy barrier E a , A is a 
pre-exponential factor, n is order of reaction and R is gas constant. If temperature rises at 
a constant heating rate, the evolution rate will first rise to a maximum value and then return 
to zero when the process of effusion ceases. The temperature (T m ) corresponding to 
maximum evolution rate determines the stability or bonding strength of the evolving 
atom/molecule. The activation energy E a for the case of constant heating rate can be 

shown to be related to T m as 

E a = -2RT m ( 3 - 2 ) 

The values of E a estimated from T m data are given in Table 3.5. Accordingly, activation 
energies for effusion of'fsiHx are - 1.05 eV and of hydrogen 1.12 and 1.28 eV. 
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3.15 Hydrogen gas evolution characteristic of PSL1 as derived from 
“peak fit” analysis of experimental data. 
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Fi^- 3.16 SiH gas evolution characteristic of (a) PSL1 and (b) PSL11 
as derived from “peak fit” analysis of experimental data. 
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Table 3.5, Thermal gas effusion data of porous silicon in the process of annealing at the 

rate of 9.6°C/min under vacuum ~ 10^ torn 


Sample 

gas evolved 

Temperature 
(T m ) for 

maximum 
partial 

pressure or 
central 
maxima (°C) 

FWHM 

Cc) 

Area under 
peak 

XI O 5 
(Torr °C) 

E a ~ -2RT m 
(eV) 

■ 

h 2 

376, 476 

38.6, 139.7 

1.2, 2.9 

1.11, 1.28 


SiH 

348 

62.8 

2.5 

1.06 


SiH 2 

345 

61.4 

3.2 

1.05 


SiH 3 

343 

66.1 

1.4 

1.05 

PSLll 

SiH 

340 

64.5 

0.4 

1.04 


SiH 2 

334 

70.4 

0.3 

1.04 


SiH 3 

333 

51.7 

0.1 

1.04 


3.5 FTIR studies 

The energies associated with the infrared region correspond to vibrational mode of 
atoms in solids. Consequently, IR spectroscopy can in principle analyze the nature of 
bonding of silicon atoms with hydrogen, oxygen, etc. Fig. 3.19 shows a typical FTIR 
spectra of PSL2 porous silicon in virgin state and after annealing upto 450°C in the 
wavenumber range of 500-3600 cm \ Also, FTIR spectra for PSL6 (virgin and annealed 
upto 250°C) and PSL13 (virgin) samples are presented in Figs 3.20 and 3.21, respectively. 
Table 3.6 summarizes the infrared absorbtion data available in the literature for vibrational 
modes and species relevant to porous silicon [5^-67 ]. The present FTIR data when 
examined carefully, it was noticed that while some peaks match well with the known 
modes (stretching or deformation)/ species, others found for the first time in virgin/annealed 
samples of porous silicon. The wavenumbers of peak positions, percentage absorption at 
peak maxima and their possible origin are listed in Tables 3.7, 3.8 and 3.9 for PSL2, PSL6 
and PSL13 samples, respectively. Also, after annealing of samples upto 250°C or 450°C, a 
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Pm- 3.19 FTIR spectra of PSL2 (a) in virgin state and 
U (b) after annealing upto 450°C. 
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Fkj 3.21 FTIR spectra of PSL 13 in virgin state. 
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Results ana Uisatssion 


1 able .> (> Summary ot the IR data available in the literature for vibrational modes and 

species/elevant to porous silicon fcp- (,y | 


s \> ■ 

I 

1 

1 

M:\ic ,-c ivies 

Characteristic peak 
maxima on band range 
wavenumber (cm' 1 ) 

Remark (s) 



814,742, 628 

Si-Si stretching 

1 

longitudinal optical and 
transverse acoustical 

mode 

(611-617)* 

" in uncorrected IR data 

2 

Si-H stretching modes of 
H* bonded hydroxyls 

3000-3800 

broad and composite 

absorption 


bending mode 

800-900 

broad and weak 

3 

I iydride species 

SiHx groups linked to Si 
stretching mode 

deformation mode 

2111 (SiH) 

2089 (SiH 2 ) 

2140 (SiH 3 ) 

628 (SiH) 

665-668 (SiH 2 ) 

906-916 (SiHj) 

862 

Wag (significant 
contribution) 

Wag scissor 
strong 

4 

Si-O-Si bonds formed 

1000-1200 

Asymmetric stretching 


with interstitial oxygen. 

(or 1106) 

broad; vibrations (appears 


Oxygen substitution into 
background of SiH 

stretching 

2143 

2240 

after ageing) 

5 

SiHx groups linked with 

2200 (0 2 -SiH 2 species) 

blue shift in PL 


oxygen: 

2258 (0 3 -SiH species) 

contribution due to Si-O-H 


stretching mode 

750 

bending mode, O-Si-O 


deformation modes 

814 

850 Oy-Si-Hx 

878 

940 

symmetrical stretching mode 
also fall in this region 

6. 

Si 02 islands 

1106 

Shift occurs to 1061 and 

asymmetrical stretching 
vibrations of bulk O-Si-O 

1200 (shoulder) 

1153 (shoulder) with ageing 


bridge 

HO-Si-H group on silica 

2240 

SiH stretching mode 


surface 

highly distorted siloxene 

2147 

include contribution from 
surface hydroxyles. 


bridges on silica 

1061 

1153 

812-814 

850 (Shoulder) 

OySiHx species and O-Si-O 

symmetrical stretching mode 
of Si0 2 layer 

7. 

Carbonate species 

1707-1720 

1288-1304 " 

zramtu 

■r h - i — 



Tabic i 7 I- HR spectrum data of PSL2 giving of peak position, percentage absorption at 

peak maxima and vv:rt>pi'iuir.g modes/species. 


\ sample 


Alter annealing upto 450 (I C 

Assignment of modes/species 

Peak position 

wavenumber 

(cm' 1 ) 

% 

Absorption 
at peak 
maxima 

Peak 

position 

wavenumber 

(cm* 1 ) 

% 

Absorption 
at peak 
maxima 

o!7 

86 

613 

84 


Si-Si stretching 

839 

i)2 




SiO-H bending mode/highly 
disordered siloxene bridge on 

- 

- 

876 

90 


deformation mode O v -Si-H 2 
species 

Broad peak 
1060 


1066 

95 


Si-O-Si asymmetric stretching 
vibration 

1527 

57 

- 

- 


New peaks; disappear after 

1560 

57 

- 



annealing 

- 

- 

1672 

70 


appears after annealing 

1720 

66 

1720 

70 


Carbonate species stretching 
mode 

21 17 

56 

- 

- 


SiH wagging 


• 

2206 

73 _ 


SiHx groups linked with 

• 

« 

2254 

17 


oxygen; C> 2 -Si-H 2 , 

** 

** 

2327 

84 __ 


Oj-Si-H species in stretching 
mode 

2860 

54 

- 

- 


New peak (disappears after 
annealing) 

3585 

75 

3585 

85 


SiO-H stretching modes of 
hydrogen bonded hydroxyls 


Table 3.8: FTIR spectrum data of PSL6 giving of peak position, percentage absorption at 
peak maxima and corresponding modes/species. 


Virein sample 


Peak position % 
wavenumber Absorption 
(cm' 1 ) at peak 


maxima 


After annealing upto 250°C | Assignment of modes/species 


% 

Absorption 
at peak 


maxima 


Si-Si stretching 
deformation mode O y -Si-Hx 
species; SiO-H bending mode 


Broad peak 


Si-O-Si asymmetric stretching 
vibration 


New peaks % absorbtion 
decreases after annealing 

Carbonate species stretching 
mode 


2144 


SiH 3 wagging; SiH stretching 
mode 


2256 


Oxygen related peaks 0 3 -Si-H 
species 


New peak; height discreases 
on annealing upto 250°C 


2858 


New peak; height disappears 
on annealing 


SiO-H stretching modes of 
hydrogen bonded hydroxyls 





Table 3.9. FTIR spectrum data of PSL13 giving of peak position, percentage absorption at 
peak maxima and corresponding modes/species. 


Virgin sample 

Assignment of modes/species 

Peak position 
wavenumber (cm' 1 ) 

% Absorption 
at peak maxima 

615 

83 

Si-Si stretching 

839 

85 

deformation mode O y -Si-Hx 
species; SiO-H bending mode 

945 

82 

SiO-H bending mode/SiH 2 scissor 

1070 

82 

Si-O-Si asymmetric stretching 
vibration 

1458 

76 

New peak 

1726 

73 

Carbonate species stretching mode 

2114 

50 

SiH wagging 

2195 

50 

0 2 -Si-H 2 stretching mode 

2247 

50 

0 3 -Si-H stretching mode 

2847 

70 

New peak 

3585 

25 

SiO-H stretching modes of 
hydrogen bonded hydroxyls 





Results and Discussion 


and some new show up. The peaks which disappear correspond to SiHx species or SiO-H 
bending mode. Such a result is consistent with the thermal gas effusion findings discussed 
in the previous section, viz., evolution of hydrogen, SiHx, etc. Also, annealing is likely to 
cause oxidation and therefore new peaks at 2254 and 2327 cm‘ l corresponding to 
stretching mode of (>2-Si-H2 and CL-Si-H species result in PSL2. It is interesting to note 
that such peaks are seen even in as prepared (i.e. virgin) PSL6. The reason for this 
perhaps lies in the longer etching time used for PSL6 and/or the peaks belong to some 
other modes/species. 

In summary, porous silicon is an inhomogenous material consisting of the pores and 
silicon skelton with nanostructures. Also, there lies a number of species on surface and/or 
in pore regions, viz., SiHx, hydrogen, hydroxyls, oxygen related species, carbonates, Si0 2 , 
etc. On annealing, hydrogen and its related species evolve and emergence of oxygen 
associated species occurs. Also, the photoluminescence (PL) is found to decrease with 
annealing. This means that non-radiative recombination centres are enhanced on the 
surface, as evident from EPR studies (see, e.g., section 3.3). The origin of PL in porous 
silicon lies with the radiative recombination centres associated with the species present on 
surface or in pore regions besides quantum confinement effects (i.e. y band gap 
enhancement) in nanostructures of silicon. Since, PS exhibits high porosity (~ 70-80%) 
and, in t rn, specific surface area of several hundreads m 2 per gram [6&], the contribution 
to PL from these regions is expected to be quite significant. This inference derives support 
even from the X-ray diffraction studies (see, e.g., section 3.2) which continue to reveal the 
presence of nanosized silicon crystallites in annealed porous silicon samples. 



Chapter 4 


Conclusions 

(1) Photoluminescence grade porous silicon has been prepared by electrochemical 
etching of p-type silicon wafers (resistivity 0.3- 1.5 fl-cm, <100> orientation) in 1:1 
mixture of 48% HF and C 2 H 5 OH at a current density of 5-10 mA/cm 2 for 15-120 minutes. 

( 2 ) Raman spectra of porous silicon contains three peaks centred around 519-521 cm' 1 , 
502-513 cm * 1 and 473-490 cm' 1 (weak and broad) originating from c-Si, nanocrystallites 
(• 5 run) and a-Si, respectively. Further, asymmetry nature of Raman spectrum towards 
lower wavenumbers suggest the existence of nanosized silicon crystallites/strain in PS. 
Annealing of PS makes the Raman spectrum progressively symmetrical as the crystallites 
gr ow in size and/or strains gradually disappear in the process. 

(3) The crystallites present in PS are of diamond-type silicon (lattice parameter 
a 5.4307 A), have nanometer size (< 5 nm) and maintain the orientation [100J of the 
original silicon wafer. Annealing of PS at 250-650 (, C leads to increase in crystallite size by 

25-50% and causes gas effusion. 

( 4 ) HPR signal with g-value 2.0061-2.0067 results due to trigonally symmetric silicon 
dangling bonds along *ci 1 1> direction. The enhancement of EPR signal in PS upon 
annealing can be attributed to (i) increase in the number of dangling bonds passivated by 
hydrogen or other species earlier and (ii) many told increase of the corresponding spin 
density. T he defects associated with the dangling bond are non-radiative type as PL 
decrease sharply with rise in the spin density. 


(5) Porous silicon effuses molecular hydrogen, oxygen, SiH x and fluride species on 
annealing in vacuum Hydrogen trapped in pores and present interstitially in silicon gets 
evolved around 376° and 476°C whereas, SiH, Sift and SiH, species effuse at -330- 
350»C; their activation energies being 1.11, 1.28eV for hydrogen and 1.05 eV for Sill, 


species. 



(b) The absorption peaks in FTIR spectrum of porous silicon either match very well 
with the known vibration modes/species, such as Si-Si stretching (-613-617 cm' 1 ), SiO-H 
stretching or bending/highly disordered siloxence bridge (~ 839 and 3585 cm' 1 ), Si-O-Si 
asymmetric stretching (- 1060-1070 cm' 1 ), carbonate species stretching mode 

( 1720-1726 cm’ 1 ). Sill* grouping (21 10-2150 cm' 1 ), or appear for the first time (e.g., at 
1458, 1527, 1570, 2330, 2840-2860 cm' 1 ). After annealing, peaks corresponding to SiHx 
species or SiO-H bending mode disappear altogether, while new ones appear at 2254 and 
2327 cm ’ 1 and correspond to stretching modes of CVSidHb and (VSi-H species 
respectively. 
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